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RESUME 

La  diffusion  d'un  faisceau  gaussien  TEMqq  par  un  objet  spherique 

est  traitee  de  fa§on  exacte  en  termes  de  fonctions  d'onde  vectorielles 

sans  aucune  restriction  de  la  dimension  ou  de  la  position  du  diffuseur. 

Les  expressions  obtenues  pour  la  puissance  absorbee  et  diffusee  sont 

presentees  sous  forme  de  combinaisons  lineaires  des  coefficients  de  Mie 

et  peuvent  8tre  des  lors  evaluees  numeriquement . Le  probleme  correspon- 

dant  a la  diffusion  d'un  faisceau  produit  par  un  laser  fonctionnant 
* 

dans  le  mode  TEMqi  est  aussi  resolu.  (NC) 


ABSTRACT 


The  scattering  of  a Gaussian  beam  TEMqq  wave  by  a spherical 
object  is  treated  exactly  in  terms  of  the  vector  wave  functions  without 
any  restriction  on  the  size  or  the  position  of  the  scatterer.  Expressions 
obtained  for  the  powers  absorbed  and  scattered  are  given  as  linear  com- 
binations of  the  well-known  Mie  coefficients  and  can  be  readily  applied 

to  numerical  computation.  The  corresponding  problem  for  the  scattering 

* 

of  a beam  produced  by  a laser  operating  in  the  TEMqj  mode  is  also 
solved.  (U) 
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1.0  INTRODUCTION 

In  infrared  system  design  and  optical  countermeasures  it  is 
important  to  have  available  quantitative  information  on  the  extinction 
properties  of  aerosols,  both  of  military  and  non-military  origin,  in  the 
atmosphere.  To  obtain  the  extinction  coefficient  and  the  scattering 
phase  function,  one  of  the  most  basic  physical  quantities  to  be  measured 
is  the  aerosol  particle  size.  For  this  purpose  light  scattering  from  a 
laser  beam  is  a technique  commonly  used.  The  standard  theoretical  tool 
in  deducing  the  size  of  a spherical  scatterer  from  the  experimental 
result  is  the  Mie  theory  developed  at  the  turn  of  the  present  century. 

To  ensure  a sufficiently  strong  scattered  signal,  the  incident  electro- 
magnetic field  is  often  in  the  form  of  a collimated  laser  beam,  whose 
characteristic  width  may  be  comparable  with  the  size  of  the  scatterer. 

In  such  cases  the  assumption  in  the  Mie  theory  [1]  that  the  incident 
field  is  an  infinite  plane  wave  is  evidently  not  valid.  The  scattering 
of  a Gaussian  beam  has  been  considered  by  Morita  et  al  [2]  but  they  have 
assumed  that  the  scatterer  was  small  compared  to  the  beam  waist.  A more 
satisfactory  approach  to  the  problem  has  been  given  by  Tsai  and 
Pogorzelski  [3].  Using  the  cylindrical  vector  solutions  of  the  wave 
equation,  they  studied  the  scattering  of  a Gaussian  beam  from  a laser 
operating  in  the  fundamental  mode  by  a spherical  particle  of  arbitrary 
size.  The  particle,  however,  is  assujned  to  be  on  the  beam  axis.  In 
many  physical  applications  - aerosol  particle  size  measurement,  for 
example  - it  is  important  to  know  the  effects  of  the  variation  in  intensity 
of  the  laser  beam  in  a given  sampling  volume  to  properly  design  the 
experiment  as  well  as  to  understand  the  measurements.  Because  of  the  lack 
of  a more  general  theory,  the  interpretation  of  the  results  of  a recent 
experiment  [4]  has  to  be  limited  to  scatterers  on  the  beam  axis. 
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In  this  report  we  present  theoretical  results  for  the  scat- 
tering of  a Gaussian  beam  produced  by  a laser  operating  in  the  fundamental 
TEMqo  mode.  Some  of  the  preliminary  results  have  been  presented  in  Ref.  5. 
The  spherical  scatterer  is  arbitrarily  situated  and  its  size  is  also 

arbitrary.  For  the  measurement  of  aerosol  particle  size  by  intracavity 

★ 

scattering  [6,7]  the  use  of  a beam  in  the  TEMgi  mode  has  been  suggested. 

The  power  scattered  and  absorbed  from  such  a beam  by  a spherical  particle 
will  also  be  given.  This  work  was  performed  at  DREV  during  1976  under 
PCN  33A11  (formerly  PCN  15B34),  Aerosols  Studies. 

2.0  EXPANSION  OF  THE  ELECTRIC  FIELD  OF  A GAUSSIAN  TEMgo  BEAM 

Before  discussing  the  more  general  theory  we  briefly  sketch  the 
results  of  Ref.  3.  It  is  assumed  that  the  intensity  of  the  electromagnetic 
beam  I(z=-Zp)  measured  at  a plane  z=-Zq  normal  to  the  beam  axis  is  given  by 

I(z=-Zq)  <r  (1) 


and  the  corresponding  electric  field  is 

^(z=-Zq)  = . (2) 

Since,  as  pointed  out  by  Carter  [8],  (2)  does  not  satisfy  the  Helmholtz 
equation  except  when  r is  small  compared  to  w^,  it  gives  only  the  near 
axis  behavior  of  ^ in  the  plane  z=-Zq.  In  fact,  Tsai  and  Pogorzelski  [3] 
have  shown  that  the  electric  field  satisfying  (2)  near  the  beam  axis 
can  be  derived  from  a Hertz  potential: 
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with  ? = -p  I-  ^ X n (4) 

In  (3)  the  cylindrical  coordinates  (r,  (Ji,  z)  refer  to  a frame  with  its 
origin  on  the  beam  axis,  h = k /l-(X/k)^  and  k^  = + jyou,  where 

a,  e and  a are  respectively  the  permittivity,  permeability  and  con- 
ductivity of  the  medium  of  propagation.  Except  that  we  have  chosen  a 
beam  wave  propagating  in  the  positive  z-direction,  the  symbols  used  here 
are  the  same  as  in  Ref.  3.  With  (3)  the  magnetic  field  can  also  be 
written  down  from 


H = V X V X n 


(5) 


In  fact,  the  electric  field  ? given  by  (3)  and  (4)  can  be  written  as  a 
linear  combination  of  the  first-order  cylindrical  wave  eigenfunctions 


Sou 


I OC 

= "-0  f 

Jo 


-x^w  2/4 

0 ’"olX  ® 


x2dX 


(6) 


To  solve  the  scattering  problem  with  the  spherical  particle  on  the  beam 
axis  Tsai  and  Pogorzelski  [3]  established  the  following  relationship 
between  and  the  spherical  wave  eigenfunctions 


Jhz 


2)1+1 

)l()l+l) 


(H-l)!  .)l-l 
C)l+1)!  ^ 


dPf  ^ (cosa) 
da 


Pi ^ (cosa) 
sina 


C7) 


The  scattered  fields  were  then  determined  through  the  continuity  condi- 
tions on  the  surface  of  the  scatterer. 
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For  the  general  problem  where  the  scatterer  is  not  restricted 
to  the  beam  axis,  we  choose  the  center  of  the  spherical  particle  as  the 
origin  O'  of  another  coordinate  system  obtained  from  the  old  system  by 
translation.  Relative  to  the  new  system,  the  origin  of  the  old  system 
is  situated  at  OCr^,  4)^ , z^)  (Fig.  1).  To  rewrite  the  Hertz  potential 
of  (3)  in  the  new  coordinates,  the  addition  theorem  for  Bessel  functions: 


Jn'f> 


(8) 


111= -CO 


can  be  used  and 


z 8a> 


cc 

I 


.-x2w  2/4  Jh(z’-2i-  = o)  X 


^ fj^^jCAr')  e^*^'  + Jg_^(Ar')  X^d\  (9) 

It  is  more  convenient  for  our  purpose  to  combine  the  cylindrical  vector 
eigenfunction  of  even  and  odd  parities  so  that 


"sX  " “osA  ■ ^ '^<•1=  n 


esA'  “sA  “osA  ^ '^esA 


Ti'e  electric  field  becomes 


t - ^'^0 


jut  r 

*^o 


.-A2a,o2/4  ^jh(z'-.z^*ZQ) 


(lO'l 


s=-® 


fm  , , + in  ...  A2d\ 

s+l,A  s-l,Aj 


(ni 
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Unlike  the  electric  field  referring  to  an  on-beam  axis  system,  the 
field  E now  contains  cylindrical  vector  waves  of  all  orders. 


To  express  (11)  in  terms  of  the  spherical  waves  n^^^ 

defined  in  a way  similar  to  (10)  it  is  necessary  to  generalize  (7).  It 
can  be  verified  that  for  any  integer  s (Appendix  A) 


jhz-^ 
e"^  m 


sX 


OD 

z 


21+1 


«,(il+l) 


(i^-s)!  .a-s+l 
(Jl+s)!  ^ 


dPj?.  (cosot) 


da 


£s 


+ —7-^ — P ^ (cosa)  n 


sina  il 


Is 


(12) 


where  P.'  are  the  Legendre  polynomials. 


Substituting  (12)  into  (11)  we  obtain 


CO  06 


= E E D^C^^.s)  m^^  + Dg(Jl,s)  n 


£=1  s=-‘» 


Jls 


(13) 


where 


D,  :j,s) 


1 A 

(Ugf^.S) 


j'V  e'^“^ 


)U)t  f 

Jn 


dXe-^V/^  X2 
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dP}  (cosa'^ 

dot 

sP;^  (co?a~ 


In  an  on-beam  axis  cdordinate  system,  E is  a linear  combination  of  the 
spherical  waves  m^^  and  n^^  only. 


3.0  SCATTERING  OF  A GAUSSIAN  TEMnn  BEAM 


To  find  the  scattered  fields  e’',  (T  produced  by  the  spherical 
particle  at  O',  the  continuity  conditions  on  the  surface  of  ’’he  scatterer 
fR  = a)  "lay  be  used 

-v  -wf 


Upj  X (E  + 2^’')  = Uj^  X 


-V  -VT  ->-t 

Uj^  X CH  + H ) = Uj^  X 


,tt 


where  S'",  H'"  arc  induced  fields  inside  the  sphere.  For  E^  a d E^  we  write 


00  00 


n 

D „ n . 

i.s  Is 


£=1 
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The  spherical  waves  m n , m and  n„  are  obtained  fm- 

XS  },S  Is  f,s 

the  corresponding  even  and  odd  parity  waves  given  in  Ref.  7.  The  coeffj 
r r t t 

cients  a b a b can  now  be  determined  from  (16)  and  (17  . 
It  can  be  shown  that  the  coefficients  are  independent  of  s such  that 


r r ,r  , r 

^ is  ^ i.s'  ’ ^ Ps  ■ ^ ^s'  ' 

^t  _ t , t , t 

2,s  ^ ?s  * * ?. s P.S  * ' ^ ' 

r r 

Furthermore,  a b arc  just  the  well-known  Mie  coeff ici '^nts  given  by 


(Appendix  B) 


r j?Cpl)  [pjg(pn  -j;(c)  [pij,(oj)1 

^ is  = • j^(pj)  [ph4hp)j'  - h4^(p)[Pij;i(Pl)] 


_ _ j ? (p)  [pi  j (.  (pi ) ] - (k'/k)ij;  (pj^)  [pj.  (d)]  

h^ihp)[Pijii(oi)]'  - (k7k)2j(i(pp[ph^i^(p)] 


where  p = ka,  = k'a  and  k'  is  the  propagation  constant  for  the  scatterer. 

Therefore,  no  scattering  coefficient  other  than  the  Mie  coefficients  in  (22' 

—^r 

and  (23)  are  involved  in  E . Dropping  the  unnecessary  subscripts  s in 


r J u r 

a . and  b , we  can  write 
is  is 


D^(ns)  . bj 
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-►  -►x* 

The  magnetic  fields  H and  H are  given  by 


OD  00 


^ i=l  S=-“ 


£S  " 


(25) 


and 


-*T 

H 


°°  ^ <”  <X) 

X.-1  S=-“>  s=-<»> 


s)  m 


(3) 

ils 


(26) 


With  the  expressions  (24)  and  (26)  for  the  scattered  fields 
we  can  write  down  at  once  a generalized  expression  for  the  power 
scattered  in  a given  direction. 

The  rate  at  which  energy  is  being  scattered  by  the  particle 
can  now  be  obtained  by  integrating  the  Poynting  vector  corresponding  to 
and  over  the  surface  of  a large  sphere  of  radius  R centered  at  O'. 
In  fact. 


lim 

W = n 

s 


f f "r  - 

-'0  Jo 


in0d0di}i 
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In  deriving  (27),  we  used  the  following  relations 


s:  r(' 


, IT  Ztt 


■"Ce  sineded^  = o 


3)*  \ 

-s-e/ 


("is*  "rs*e  - "Ise 


f.fj 


' (3)  (3)‘ 

"“use  "ji's' 


(3)*  ^ 


in9d6d({! 


- K X 

■ 1,1'  s,s'  21*1 


In  a similar  manner,  the  rate  at  which  energy  is  being  scattered  and 
absorbed  by  the  particle  can  be  found  by  evaluating  the  integral 


lim 


TT  ^2ir 


R^sineCEjHj*.  eJhJ  - E,hJ'  - EjH^ldOd*  (31) 


to  give 


♦ bJ|Dg(a,s)|2j 


UNCLASSIFIED 

11 


For  the  special  case  when  the  beam  axis  passes  through  the 
center  of  the  spherical  scatterer  r^  = 0,  (14)  becomes 


^bC^s) 


+6 

s 


i)  X 


A()l,s;X) 

X^dX 

B(£,s;X) 


(33) 


After  some  straightforward  algebraic  manipulation  we  arrive  at  the  following 
results  for  on-beam  axis  scattering 


W 

s 


CD 

2(2t+l) 

ku>  / j [t(t+l)]2 
£=1 


[' 


(34) 


and 


W 

t 


2(2£+l) 

TT(£7IT]2 


|d^(a)i2 


(35) 


where 


dA(t) 


4 

+I 

•'0 


^ X^iX  (36) 


and 


dB^) 


OD 

= !0  f gjh(z^+ZQ)  Pj^cosa)  ^2^^ 

8 sina 


(37) 
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In  the  infinite  plane  wave  limit  corresponding  to  both  and 

d_(il)  approach  to  a common  value  )l(Jl+l)/2  and  the  well-known  results 
o 


and 


00 


(38) 


(39) 


are  obtained. 


4.0  EXPANSION  OF  THE  ELECTRIC  FIELD  OF  A TEMqi  BEAM 


The  first-order  oscillation  mode  in  a laser  cavity  with  cylin- 

* 

drical  symmetry  about  the  beam  axis  is  often  referred  to  as  the  TEMqi 
mode.  Its  intensity  distribution  in  a plane  normal  to  the  beam  axis  is 

_ -2r^ /w  ^ 

I(z=^Zq)  r^  e 0 (40) 

and  the  electric  field  near  the  beam  axis  is 

^(z=-Zo)  " “x 

In  contradistinction  to  the  fundamental  mode  discussed  above,  such  a field 
cannot  be  derived  from  a single  Hertz  potential.  Two  complementary  Hertz 
potentials  and  of  the  electric  and  magnetic  types  [9]  are 

required  and  the  corresponding  partial  fields  are  given  by 

It  * « 


^ x ^ X 


= c 


(42) 
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and 

^ _y  3 V i^C2)^  j5(2)  = ^ X 7 X (43) 

o c 

From  the  general  theory  of  vector  potentials,  it  is  sufficient  to  consider 
of  the  form 

= xu^  (44) 

An  immediate  consequence  of  using  a Hertz  potential  of  the  electric  type 
is  that  is  not  completely  transverse.  However,  the  longitudinal  compo- 
nent of  ^ does  not  contribute  to  I(z=-Zp)  and  it  will  be  seen  that  the 
field  ^ is  dominated  by  its  transverse  component  near  the  beam  axis. 

The  Hertz  potentials  t(),  x can  be  expressed  as 

^ cos^^,  f(z)  (45) 

X = ^ e"^“^  r^sin^i)!  e ^ ^'^0  g(z)  (46) 

The  functions  f(z),  g(z)  are  chosen  such  that  rp  and  x satisfy  the 

df  I 

scalar  wave  equation  as  well  as  the  boundary  conditions  2-  z “ 
g(z=-ZQ)  = 1.  As  in  the  case  of  the  fundamental  mode,  and  x ^ 
can  be  expressed  in  terms  of  the  solutions  of  the  scalar  wave  equation 
in  cylindrical  coordinates.  Let  us  define  the  functions  q2(X),  qQ(^) 
and  q 2(^)  to  be 


r 
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where  i=l,  2,  3 and  5^=2,  S2=0,  s^=-2.  Using  the  boundary  conditions 
of  Ip  and  X and  the  Fourier-Bessel  theorem,  we  can  obtain 

3 

4'  = ^ e-j“t+js^ij>  |j^  J — Jg^(Xr)  XdX  I 


and 


X = 


3 00 

g qsiW  J,.(xr) 


XdX 

(49) 


where  <52=52=1,  <52=2. 


The  partial  fields  and  are  then  given  by 


^(1)  ^ g-jut  js  * 6 


i=l  Jo 


AX)  X 


3Js  (Xr)  ^ s.  ^ 

^ 3r  “r  ■ Jsi(Xr)  ^ Jsi(^r)  J 


X2 


XdX  (50) 


3 * 


^ .(Xr)  ^ 

^ Jsi(Xr)  I XdX 


(51) 
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The  partial  field  i3 

transverse  but  ^(1) 

has  a longitudinal  component , 

It  is  clear  from  (50)  that  near  the  beam  axis  (r-K))  the  transverse  compo- 
nent dominates.  Both  Jd) 

and  can  now  be  expressed  in  terms  of  the 

cylindrical  eigenfunctions  so  that 


t = e-^“^ 


3 « 

z'-ff 

i=l  Jo 


h "siX  jms-x  [ . 


(52) 


5.0  SCATTERING  OF  A TEMqi  BEAM 

The  electric  field  in  (52)  is  expressed  in  terms  of  the 
on-beam  axis  coordinates.  To  find  the  fields  scattered  by  an  arbitrarily 
situated  sphere  we  again  introduce  an  off-axis  system  with  coordinates 
(r',  (j)',  z')  as  before.  By  applying  the  addition  theorem  for  Bessel 
functions,  the  field  ? in  the  off-axis  system  becomes 
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To  transform  the  cylindrical  wave  functions  m n , to  the 

^ sA  sX 

spherical  functions,  we  need  an  identity  for  n , similar  to  that  for 
^ sX 

given  in  (12).  It  can  be  shown  that  (Appendix  A) 


Jhz  ^ 21*  1 ()l-s)!  ,t-s^l,  rdP%. 

sX  ^ «,(jl+  1)  (i+s) ! •’  ^ [ 


(cosa)  ■+■ 

T HI, 

da  Is 


4=1 


" -rk- 


(55) 


and  hence 


^ S [j(^.s)  m,^^  + K(4,s)  n J 


(56) 


4=1  s=-<' 


where 


J(il.s)  = 


-jut 


3 


/ C(s-Si,i,X)  B(4.s;a)-(-1)^  A(4,s;X)]  dX 

i=l  Jo  ^ 

(57) 


and  K(4,s)  is  obtained  from  J(4,s)  by  interchanging  A(4,s;X)  with 
B(4,s;X).  Because  t in  (56)  has  exactly  the  same  f)rm  as  (13),  the 
scattered  fields  can  be  written  down  at  once 


and 


^ p4  ”^^4s  ^ ^"^4  "^^4sl 

s=-« 

S 2 p4  ^ n^^jJ  (59) 

4=1  S=-oo 
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The  rates  at  which  energy  is  being  scattered  away  from  the  beam  are 


W = 


s ~ liu 


Za(*.+  1)  (A+S)! 

2t+l  (a-s)! 


11=1  S=-“> 


jla’'^  J()l,s)l^  + K(£.s)|^j 


(60) 


and  W^,  the  rate  at  which  energy  is  scattered  and  absorbed,  is 


*t  = i S 

1=1  S=-“> 


a{t+l)  (H+s)! 
211+1  («.-s) ! 


a’'Jj()t.s)|^  + b’^jKCi.s)!^ 


(61) 


The  coefficients  b’^^  in  (58)  - (61)  are  again  the  Mie  coefficients 
given  in  (22)  and  (23). 


6.0  CONCLUSIONS 

The  results  of  Tsai  and  Pogorzelski  [3]  derived  for  a Gaussian 
beam  corresponding  to  the  fundamental  mode  TEMgo  have  been  generalized. 
In  the  first  place,  the  incident  fields  are  expressed  in  terms  of  the 
cylindrical  and  spherical  vector  wave  eigenfunctions  respectively  in 
an  off-beam-axis  coordinate  system.  This  enables  us  to  obtain  the 
scattered  fields  as  well  as  the  powers  scattered  and  absorbed  from  the 
beam  by  an  arbitrarily  situated  spherical  object.  We  have  shown  that 
although  the  incident  and  the  scattered  fields  in  an  off-beam-axis 
system  are  linear  combinations  of  spherical  wave  functions  of  all 
orders  (i.e.  m^^^,  n^^^  with  s equal  to  any  integer),  the  scattering 
coefficients  are  independent  of  s and  are  just  the  Mie  coefficients. 
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Parallel  results  are  also  obtained  for  the  cylindrical  symmetric  mode 

it 

TEMoi.  These  theoretical  results  together  with  a standard  program  for 
calculating  Mie  coefficients  [10]  can  be  readily  applied  to  nunerical 
computations  for  the  scattering  of  beam  waves  without  imposing 
restrictions  on  the  size  or  position  of  the  scatterer.  With  the  more 
general  theory,  aerosol  particle  sizing  can  be  more  efficiently  carried 
out  using  the  laser  scattering  technique. 
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we  can  show  after  rather  lengthy  but  elementary  manipulations  that 
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is  a vector  which  has  no  component  along  Uj^.  In  a similar  manner,  the 
u component  of  N can  be  proved  to  be  given  by 
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From  the  fact  that  m™  contains  no  Uq  component,  it  follows  that 
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may  only 
that  ^ 
free  (as 
a linear 


have  a non-zero  component  along  u . To  complete  the  proof 
is  indeed  a null  vector,  we  note  that  since  must  be  divergence 
fj,  n^^,  are  all  divergence  free)  it  can  only  be  expressed  as 
combination  of  and  n^^  i.e. 
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In  (A-9)  only  n^^  has  Uj^  component  and  hence  6^=0  following  from 

the  linear  independence  of  the  vectors  Similarly,  since  m^^  has 

■*  -V'  ' 

components  while  N does  not,  the  coefficients  must  also  vanish. 
This  proves  (A-2)  and  (A-1)  can  be  obtained  at  once  from  (A-2)  and  the 
fact 
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APPENDIX  B 


'^R  ^ ^Is  = ^ Hs  = ^^s 


From  the  boundary  conditions  (16,  17)  on  the  spherical  surface  R - a, 
we  obtain  the  following  equations: 


l:  L (s^  * RW-  |1  a‘^ 
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In  eqs  (B-10)  and 


0 s ^ 0 

, for  „ 

1 s = 0 


since  the  two  sets  of  vectors  {^£s’  ^Is^'  ^Is^} 

have  respectively  the  same  angular  dependence.  Relations  similar  to 
eqs  (B-9)  - (B-11)  can  readily  be  written  down  for  any  two  of  the  six 
vectors.  From  these  relations  and  eqs  (B-1)  and  (B-2) , we  obtain  that 

j£(ka)  = a^gh^^^Cka)  - a^^j.t(k'a)  = 0 

[pj£(p)]  + a^g[ph|^^ (p)]  - a^g[Pjj£(Pj)]  =0 


j£(p)  + bjgh^(p)  - bjg  j£(p^)  = 0 


where  p = ka,  Pj^  = k'a. 
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It  is  clear  that  the  unknowns  a»  , b 
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equations  are  independent  of  the  values  of  s.  Furthermore,  for  any 
arbitrary  integers 
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[pj£(p)]'  + bj^  {Phf^(p)]'  - ^b*g  [Plj^CPj)]’  = 0 (B-15) 
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